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Abstract 

The seesaw mechanism is conceived on the basis that a mass scale, £, and a dimensionless 
scale, s, can be fine-tuned in order to control the dynamics of active and sterile neutrinos 
through cosmon-type equations of motion: the seesaw cosmon equations. This allows for 
sterile neutrinos to be a dark matter candidate. In this scenario, the dynamical masses 
and energy densities of active and sterile neutrinos can be consistently embedded into the 
generalized Chaplygin gas (GCG), the unified dark sector model. In addition, dark matter 
adiabatically coupled to dark energy allows for a natural decoupling of the (active) mass 
varying neutrino (MaVaN) component from the dark sector. Thus MaVaN's turn into a 
secondary effect. Through the scale parameters, £ and s, the proposed scenario allows for a 
convergence among three distinct frameworks: the cosmon scenario, the seesaw mechanism 
for mass generation and the GCG model. It is found that the equation of state of the 
perturbations is the very one of the GCG background cosmology so that all the results from 
this approach are maintained, being smoothly modified by active neutrinos. Constrained 
by the seesaw relations, it is shown that the mass varying mechanism is responsible for the 
stability against linear perturbations and is indirectly related to the late time cosmological 
acceleration. 

PACS numbers: 98.80.-k, 12.10.-g, 14.60.St, 95.36. +x 



*Electronic address: alexeb@ufscar.br, alexeb@ifi.unicamp.br 

'''Also at Instituto Superior Tecnico, Departamento de Fisica, Av. Rovisco Pais, 1, 1049-001, Lisboa, Portu- 
gal, as a visiting scholar 
■^Electronic address: orfeu@cosmos.ist.utl.pt 



3 Also at Instituto de Plasmas e Fusao Nuclear, 1ST, Lisbon 



1 



I. INTRODUCTION 



Theoretical efforts to interpret the observational data and to understand the nature of the 
dark sector do necessarily involve a vivid interplay between general relativity, astrophysics 
and particle physics. Since the simplest solution to account for the late time accelerated 
expansion of the Universe, the one given in terms of a tiny positive cosmological constant, 
is plagued with conceptual problems, one has been compelled to examine different solutions 



ft 



Motivated by the high energy physics, an interesting alternative for obtaining the neces- 
sary negative pressure to account for the accelerated expansion involves the dynamics of a 
scalar field, 0, evolving slowly down its potential V(4>) js, 9|. These models assume that the 



vacuum u v can vary; a feature discussed muci earlier (lOj. O ;her alternatives include 



fc-essence 



11 



121 ]. phantom energy models [13|, Il4j . 



cosmon 



fields [15|, |l6| , and also several 



17H19|. 



modifications of gravity 

A challenging related issue concerns models of mass varying particles j2ol[22|. These are 
coupled to a light scalar field that drives their energy through their dynamical mass. The 
idea of this mass varying mechanism [?], 22 1 is to introduce a coupling between a relic 
particle, usually neutrinos, and the dark sector: dark energy or dark matter separately, or all 



together 



23 



24j . Such models admit an adiabatic regime in which the scalar field, usually 



related with dark energy, rolls down the minimum of its effective potential, which is set 
by the dark matter dynamical mass. As a direct consequence of this new interaction, the 
particle mass is altered by the dynamics of the scalar field. 

Due to phenomenological reasons, one still expects a small contribution from neutrinos 
to the cosmic dynamics. In fact, it is well-known that the active neutrino masses are tiny 
as compared to the masses of the charged fermions. The smallness of the neutrino masses is 
usually understood in terms of the seesaw mechanism in extensions of the standard model 
(SM) of the electroweak (EW) interactions. The EW interactions involve only left handed 
neutrinos such that no renormalizable mass term for them is compatible with the SM gauge 
symmetry SU(2)l £g> U{l)y. Once one assumes the conservation of baryon and lepton num- 
bers, the seesaw mechanism admits neutrino masses from dimension five operators. Neutrino 
masses should then involve two powers of the vacuum expectation value of the Higgs dou- 
blet. These masses are suppressed by the inverse power of a large mass scale, M, of a right 
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handed Majorana neutrino. This super-heavy Majorana neutrino is associated to lepton 
number violating effects in extensions of the SM. 

Assuming that the sterile neutrino mass, M. = M.{4>), exhibits a dynamical behavior 
driven by the scalar field, 0, then the sterile neutrino becomes an interesting candidate 
for the aforementioned mass varying dark matter. This, through the seesaw mechanism, 
gives origin to a mass varying neutrino (MaVaN) 2CH22] component corresponding to active 
neutrinos that exhibit secondary mass effects due to the seesaw mechanism. In models of 
coupled dark energy, in which the scalar field couples to other matter components, it is 
natural to expect a coupling to active and to sterile neutrinos. 

In this work, one explores the consequences of such a coupling. This leads to cosmological 
scaling solutions where dark energy tracks the evolution of matter and/or radiation. For 
the present cosmological epoch, it is predicted that the energy densities of dark energy and 
matter are of the same order of magnitude. If the tracking can be altered by the growing 
mass of neutrinos, such that they become nonrelativistic (NR) at low redshift, one can 
interestingly match the observed universe. 

One should notice that mass varying dark matter is unusual in the formulation of the 
MaVaN models. The unification of dark energy and dark matter naturally offers this pos- 
sibility. The generalized Chaplygin gas (GCG) is particularly relevant in this respect [2|-[4| 



as it is shown to be consistent with the observational constraints from CMB 



25], supernova 



26H28], gravitational lensing surveys 29], and gamma ray bursts 301 ] . 

Once one sets the mass dependence on 0, the dark sector scalar field, one notices that 



the cosmological evolution of the unified fluid composed by mass varying dar 



i matter and 



evolving dark energy has a dynamics similar to that of the cosmon field [16|, |31|, [32] . Thus, at 
least partially, it turns out that the mass varying mechanism is the main agent of stability 
and of the cosmic acceleration. Any cosmological fluid which effectively reproduces the 
effects of the GCG naturally offers analogous possibilities 33]. 

Assuming that the scalar field drives the behavior of the masses of active and sterile 
neutrinos, the mass varying mechanism and the conditions for a stable cosmological scenario 
naturally emerges in the context of the GCC model. Remarkably, the seesaw masses and 
energy densities of active and sterile neutrinos can be consistently embedded in the GCG 
scenario without any additional assumption. Such a connection is mediated by a mass scale, 
£, and a dimensionless scale, s, similarly to cosmon-type dynamical equations 16|, |3l|, 132 ]. 
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These are dubbed as seesaw cosmon equations. This procedure gives origin to a remarkable 
convergence of three distinct frameworks: the cosmon-like dynamics, the seesaw mechanism 
for mass generation and the GCG scenario. The equation of state of the perturbations is 
the very one of the background cosmology so that all effective results arising from the GCG 
model are maintained, although modified by neutrinos in a quite subtle way. 

This work is organized as follows. The decoupling mechanism for active and sterile neu- 
trinos is described in section II. This is obtained through the coupling of a scalar field in the 
seesaw relations, which leads to the dynamical properties ensued by the mass varying mecha- 
nism. In section III, the interplay with the GCG model is discussed. In section IV, the main 
properties of a unified treatment of dark matter and dark energy in a dark energy scenario 
for a A-like equation of state, p = — p, is reviewed. This leads to a dynamical mass prescrip- 
tion different from the one obtained through the simplified version of the seesaw mechanism 
explored in section III. Specifical mass dependencies on the scalar field for which the mass 
varying dark matter results in an effective GCG model are discussed. Finally, the results 
for energy densities, growing neutrino mass and stability conditions from a positive squared 
speed of sound c\ are discussed in section V. Three cosmological scenarios are examined, actu- 
ally associated to three different growing mass relations: p(4>) = 0exp [—3(a + 1)0/2] (Case 
01), /i(0) = exp [-3(a + 1)0/2] (Case 02) and /i(</>) = tanh [3(a + 1)0/2]^ exp [-3(a + 1)0] 
(Case 03). The obtained results indicate that the proposed approach is quite appealing as 
it unifies neutrinos, dark matter and dark energy. In section VI, the main implications of 
the developed scenario are summarized. 



II. THE SEESAW MECHANISM FOR MAYAN'S 



The connection of mass varying dark matter with neutrinos does provided, as will be 
seen, interesting constraints on the neutrino masses, on the dark energy density, and on 
the equations of state and the stability conditions. This can be understand through the 
equations arising from the Lagrangian densities of active (A) and sterile (S) neutrinos, ip AjS , 

C A ,S = i$A,8lv.& 1 ll>A,3 + k A ,si>A,S^A,S, (1) 
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where two mass scales, k A = p and k s = A4, have been introduced. The seesaw mechanism 
suggests that p is small, while M. should be large: 



- p = (M/2)[l - y/l + 4(m/M) 2 ] and M = (M/2)[l + y/l + 4(m/M) 2 ]. (2) 
These states correspond to the eigenvalues of the mass matrix 



written in the 



m 
m M 

orthogonal basis of chiral left- and right-handed neutrinos, z/£ , related with the matter 
fields, ip AjS , by 

= Vs 2 + 1(sz/l - v R ) and ip s = \/ s 2 + + sv K ), (3) 

where the dimensionless quantity s = \fpjM. has been introduced. It follows that C A +C S = 
£j j + £ h + £j jH , where £^ correspond to the kinetic terms and £j Ji? yields the mass mixing 
terms. The equivalence between the stress-energy tensor in the chiral basis and the matter 
field basis follows from the relationships for the energy density p A + p s = + p R and for 



the pressure p A + p s = + p R . After introducing an auxiliary mass scale £ = y/pM. = m, 
one can define two energy scales, p 1 = (p s + p A )/2 and p 2 = (p s — p A )/2 which evolve as 
reciprocally coupled cosmon-type equations. In a FRW universe it corresponds to 

Pi+3 „ (Pl+p0 -^ + ^0, (4, 

and 

p 2 + 3F(p 2+ p 2 )-^^ + 0^^ = O, (5) 

where H = a/a is the expansion rate of the universe and the overdot denotes differentiation 
with respect to time (' = d/dt). The third terms in the above Eqs. correspond to a 
mass varying term (see Eq. f l32]) in the Appendix), while the last ones are associated to the 
exchange of energy due to the non-adiabatical behavior of energy densities 1 and 2. Notice 
that there is no sense in defining a Lagrangian density for the component p 2 , which does not 
correspond to an observable energy density scale. The p x component can be identified to 
matter fields, while p 2 is an auxiliary energy density which measures the coupling between 
the physical observables. In fact, both energy densities, 1 and 2, are driven by cosmon-type 
equations. The cosmon signature is revealed by the dependence of the scales £ and s on 



the value of a slowly varying classical scalar field 4>, the seesaw cosmon field [16| . Since the 
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scalar field depends on the scale factor a, </> = 0(a), the seesaw mass terms get transformed 
into dynamical quantities, p(4>) and M.(4>). After combining Eqs. ([2]), and observing that, 

p A {a,£,,s) = p A (a,£s), (6) 

and 

Ps(a,Z,s) = p s (a,Z/s), (7) 

it is easy to identify the evolution of these with the evolution of active and sterile energy 
densities, p A and p s , by means of decoupled equations 

p A + 3H(p A+ p A )-^^ = 0, (8) 



and 



as 



d0 w ' V d0 

and 



A + Sff(» + »)-#^ = 0, (9) 



dlnf . /dlnu dln.M 
(V2) + — TT- 



dins /Mln/i dlnA^ 



Thus, mass varying mechanism translates the dependence of the mass terms on the scale 
factor, a, i.e. p(a) and /A (a) (see the Appendix). The coupling between relic particles 
and the scalar field as described by Eqs. (JS}-© is relevant only for NR fluids . Since the 
strength of the coupling is suppressed by the relativistic pressure, as long as particles are 
ultra-relativistic (UR), the active and sterile neutrino energy densities, p A and p s , decouple 
from each other and evolve adiabatically (35], remaining coupled only to the scalar field. 
That is, in the UR regimes, J£ ~ |f oc (p - 3p) « 0. 

In a previous work, it has been suggested that one could treat MaVaN's as a perturbative 



component derived from an unperturbed adiabatic energy density solution p^ [23|, |24( . The 
above results provides a quantitative justification for that. As can be seen, all the information 
from the dark sector (dark energy plus dark matter) acting on the (active) neutrino sector is 
carried out by the explicit dependence of p = p{4>). From the cosmological point of view, it 
results in mass eigenstantes for active award sterile neutrinos that evolve separately, which is 
not the case for the coupled chiral eigenstates z/£ and v n . At primordial times, when s 2 ~ 1, 
such mass eigenstates are indistinguishable, and the chiral eigenstates are well-defined. At 



late times they turn into sterile and active mass eigenstates, maintaining the identity of the 
flavour sectors. 

In order to proceed, one identifies the large mass energy density component of the above 
equations, p s , to the energy density of dark matter, then the seesaw cosmon framework 
provides the connection between dark matter and dark energy through the cosmon field 
equation l]], Q , 

H + 3H(p, + P^) + ^^= 0, (10) 

originally written as 

with the usual assignments for the p^ and (c. f. Eqs. flTTl) ). Eqs. - (fit)]) result in the 
adiabatic equation for the dark sector, 

p DS + 3H(p DS +p DS ) = 0, (12) 

with p DS = p<f> + p s = p<j> + p x + p 2 and H 2 = p DS (with H in units of H and p DS in units 
of Pent = 3H 2 /8nG). Despite the intrinsic dependence on (ft, the equation of motion for 
the dark sector is not modified by p A , the active neutrino energy density component. The 
cosmological dependence of p A on (f> can be computed through Eq. (jSJ) considering active 
neutrinos test fluid. 

The phenomenological consistency of the proposed scenario can be assessed quantitatively 
expressing p A and p s as energy densities of a degenerate fermion gas (DFG) at different 
relativistic regimes, 

p A (a) = (8^ 2 )-V(«) 4 [77(a) (277(a) 2 + l)vVa) 2 + 1 - arc sinh (77(a)) 



p s (a) = (87r 2 )-V(a) 4 ^ 8 7S 2 77(a)(27 2 s 4 77(a) 2 + 1) y/-y*s*rj(a) 2 + 1 - arcsinh (js 2 7]{a)) , (13) 

where 77(a) = (T q A )/(a p(a)) and the relation between the fluid thermodynamic regimes 
is parameterized by the coefficient 7 = q s /q A . As discussed in the Appendix, the DFG 
prescription is suitable for describing the transition between UR and NR regimes. In this 
case, the effects due to the coupled dark matter {M.(<t>)) and dark energy ((f)) can be monitored 
through Eqs. ©-flT^. 

In the Fig. [TJ one can see the exact correspondence between the abovementioned energy 
densities and the "modified" scale parameter 7s 2 . In the NR limit of a DFG, one has 
Pa/Ps ~ T 3 s 2 . 



The characteristic magnitude of the active neutrino masses involves an appropriate combi- 
nation of dimensionless Yukawa couplings, Y d , p 3 ~ m 2 /M with m 3 - ~ Y$ v (v ~ 2 x 10 11 eV). 
Consistency with the observed flavour oscillations requires for the neutrino mass at least 
one neutrino to have fa > 0.05 eV. That is, Y j of the order one implies an upper bound 
for the large mass scale M. < 10 23 eV, from which follows that s 2 > 10~ 24 . Given that the 
present value of the ratio p A j ' p s is, for phenomenological reasons, of O(10~ 2 ), a DFG of 
active neutrinos, at least approximately in the NR regime, leads to 7" 3 ~ 10 22 , and hence, 
two widely different momentum scales for active and sterile neutrinos, q s /<iA ~ 10" 7 . Thus, 
sterile neutrinos behave like ultra cold dark matter (CDM). In addition, if one assumes 
that p^ could be parameterized by a A0 4 theory, or some type of quintessence potential, the 
adiabatic evolution of the scalar field allows for assigning a mass for the scalar field. Such 
predictions for are consistent with the lower bound on the mass derived from assump- 
tions about the adiabatic evolution of the scalar field since nucleosynthesis. The scalar field 
mass should then be greater than the rate of expansion at nucleosynthesis, which is of order 
of ~ 10 -16 eV. However, the conditions for the adiabatic regime for the light scalar field to 
be settled down at the minimum of its potential prior to nucleosynthesis are quite model 
dependent. 

The above discussion shows that a considerable fine-tuning is required to generate the 
different scalar field mass, a issue that demands for an embedding of the model in a more 
encompassing framework, such as, for instance, the Minimal Supersymmetric SM with the 
addition of one singlet chiral superfield 36]. Another possible avenue involves a unified 
picture of dark energy and dark matter where the former corresponds to an additional 



singlet scalar field which one identifies to the quintessence field 34|, |37|, while the latter 



with the quantum excitations of this singlet, which can be coupled to the Higgs boson. 



III. THE INTERPLAY WITH THE GCG MODEL 



One considers now some generic features of the proposed model and explores its connec- 
tion with the GCG model. 



The GCG model is characterized by an exotic equation of state 4], |25[] given by 



p = -A s p 



Po 



(14) 



S 



where A s and a are constants. This equation of state can be obtained from a generalized 
Born-Infeld action |4|. Several studies yield convincing evidence that the GCG scenario is a 
phenomenologically consistent approach to explain the accelerated expansion of the universe. 
Inserting the above equation of state into the unperturbed energy conservation Eq. (|12|) . one 
obtains through a straightforward integration [4j 



P = Po 



A s + 



:i-a. 



and 



P 



-A s p 



A, 



(a/a ) 3(l+a) 

(I -A.) 



l/(l + a) 



-a/(l+a) 



(15) 



(16) 



(a/a ) v \ 

One of the most striking features of the above equations is that the energy density interpo- 
lates between a dust dominated phase in the past, where p oc a~ 3 , and a de-Sitter phase, 
p = —p, at late times. This evolution is ruled by the model parameters, a and A s , which 
are positive and smaller than unity. Of course, a = corresponds to the ACDM model. 
It sets < a < 1. For A s = 0, GCG behaves like matter, whereas for A s = 1, it behaves 
as a cosmological constant. Hence, in order to consider it as unified model for dark matter 
and dark energy one has to exclude these two possibilities so that A s must lie in the range 
< A s < 1. This property makes the GCG model an interesting candidate for the unifi- 
cation of dark matter and dark energy, i. e. for the dark sector energy density p DS of our 
proposal. 

Th e GCG can be described by an underlying scalar field, <ft, which can be either real 
or complex 0, 4]. In the former case, one can identify it with the cosmon field <p so 
that p and p are given in terms of by: 

1 •„ 



P 
P 



(17) 



This allows for obtaining the dependence of the scalar field, 0, on the scale factor, a, and the 
explicit dependence of p, p and V on cf). Following Ref. 27J, one obtains through Eqs. (flB"]) - 

m 



2 («) 



Pop. - A.) 

{a/a 



|3(«+l) 



A s + 



A, 



- -a/(a + l) 



(a/ ao) 



3(a+l) 



(18) 



and assuming a flat universe described by Friedmann equation H 2 = p (again with H in 
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units of H and p in units of p Crit = 3H 2 /8tiG), one gets 

1 



0(a) 



3(a + l) 



In 



i-A.(i-(o/ao; 



3(a+l)\ 



yjl - A.(l - («/«o) 3(q+1) ) + V^Aj ' 



(19) 



where it is assumed that 



(a = ao) 



1 



3(a + l) 

One then readily finds the scalar field potential, 



In 



(20) 



V(<f>) = ^ Ai+ Q p { [cosh (3(a + 1)0/2)] + [cosh (3(a + 1)0/2)] ^ 



(21) 



Thus, the dynamics of the GCG scalar field is given by the evolution of on the above 
potential. 

Turning now to the simplest version of the seesaw mechanism for which a single flavour 
neutrino mass is linearly related to the scalar field, 0, i. e. m(4>) ~ 0, after observing that 
the logarithm of the squared scale parameter s 2 has an analytical structure similar to that 

of 0(a), 

~y/l + 4(m/M) 2 - 1 



ln(s 2 ) = \n(p/M) = In 



(22) 



y/l + 4(m/M) 2 + 1 

one can rewrite the auxiliary scale m/M in terms of the GCG parameters A s and a, and on 
terms of the scale factor, a, as 



m(a) m a(a + i) l 
-a 2 



.4, 



[a/a ) 



i(* + i) 



(23) 



M(a) M 2 V 1 - A s 

where a is determined by a phenomenological adjustment (it can be set equal to unity if 
one redefines A s and p in Eq. ( I15p ). After a simple mathematical manipulation, one obtains 



In (s 2 ) = In (p/M) = -3(a + l)0(a), 



Eqs. ( 122]) - (124[) yield the following relationships for active and sterile neutrino masses 
p{(j)) = 0exp[-3(a + 1)0/2] and M{4>) = 0exp [+3(a + 1)0/2]. 



(24) 



(25) 



The behavior set by Eqs. (1231) and (1251) implies to the scale parameter s an exponential de- 
pendence on 0, s = exp [— 3(q + 1)0/2], corresponding to an exponential divergency between 



mass scales. This is consistent with several classes of quintessence models 



39|. Naturally, 
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once the prescription for masses and couplings is known, the behavior of the neutrinos can 
be understood through Eqs. (jHJ) and (jUJ). 

It is evident that the procedure above described is fairly general. The form of M{4>) and 
modifications on the equation of state for dark energy can lead to quite different scenarios. 
For instance, it also admits constant mass dark matter, which should result in neutrino 
masses given by //(a) oc exp [— 3(a + 1) 0(a)], or even a constant mass scale £ = m, as will be 
discussed in the analysis of the following section. 

IV. MASS VARYING DARK MATTER FROM AN EFFECTIVE GCG 

Aiming to obtain a deeper understanding of the proposal of the previous section, one 
shows that the coupling of dark matter to dark energy through the mass varying mechanism 
in the seesaw cosmon scenario can be indeed matched with the GCG model. In order 
to perform that, one assumes that the GCG equation of state Eq. ( TH|) just describes the 
behavior of a fluid with energy density p DS . 

From Eqs. ([9]) and ffTUj) . one sees that if the seesaw cosmon field obeys the following 
equation of state 



thus the problem consists in obtaining the relationship between the scalar potential U{4>) 
and the variable mass Ai(4>) which satisfies the equation 



a stationary condition directly obtained from Eq. ( TTU1) . The effective potential U EB (4>) that 
sets the evolution of the scalar field has two terms: the first one arising from the original 
quintessence potential U(<f>), and the second one due to a coupling to the dynamical mass 
M.((f>). For suitable choices of U{4>) and for couplings satisfying Eq. f )27|) . the competition 
between these terms leads to a minimum of the effective potential. In the quasi-staXic 
regime, it is possible for the field to adjust itself to the minimum of the potential in an 
adiabatic way. The timescale for that must be short when compared to the timescale over 
which the background density is changing. In this regime, matter and scalar field are tightly 
coupled and evolve effectively as a single fluid. In the proposed approach, once assuming 




(26) 



dj/ffl , dps dM(4>) = dU Eft (^) 
d0 dM d0 d0 







(27) 
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the equation of state p^ = —p</>, the described behavior is rather natural. Eq. ( 126|) leads to 
Pds+Pds = Ps+Ps which, in the CDM limit, gives p s = and p DS +p D s = Ps(p) = M(a)n(a), 
where n(a) is the particle number density. Since the dependence of Ai on a is mediated by 
0(a), i. e. Ai(a) = M((j)(a)), from Eqs. ffTBT) . ([To]) and f[T§j) . it follows after some mathematical 
manipulations that 

2a 
1 + a 

(28) 

which is consistent with Eq. f l27|) . Hence, one sees that the adequacy to the adiabatic regime 
is conditioned by the mass varying behavior. 

In a previous study |38|, it has shown how the behavior of ftA{4>) and U(4>) differs from 
the one of 0(a) and U{4>) in the GCG model, and how the composed fluid deviates from 
the GCG scenario. It has been assumed that the mass varying dark matter behaves like a 
DFG in the relativistic regime (hot dark matter (HDM)) and in the NR regime (CDM). For 
the mass varying CDM coupled to dark energy with p^ = — p^, the GCG leads to similar 
predictions for w = p^/p^, independently of the scale parameter a. The same is not true for 
HDM as the GCG- like behavior holds just close to the present (a ~ a ). 

V. ENERGY DENSITY, NEUTRINO MASS AND STABILITY 

As mentioned, the evolution of matter components are constrained by several factors. The 
precise peak location of the CMB anisotropies, the change in the growth of cosmic structures, 
and the properties of nonlinear structures which provide detailed information to test such 
models. In here a preliminary analysis is performed to compare neutrino masses and the 
corresponding conditions for stability for three cases of neutrino dynamical masses. Due to 
previous arguments, these cases take place in the context of the GCG model. The first two 
cases correspond to a neutrino mass dependence on the cosmon field. For the first case, the 
off-diagonal mass matrix elements, the Dirac mass terms, Tn(4>), have a linear dependence 
on 0, m{4>) ~ 0, which results in a tiny mass eigenvalue p{4>) = 0exp [— 3(q + 1)0/2]. For 
the second case, the off-diagonal mass matrix elements are constants, m(0) ~ const., and 
the dynamical masses are generated by the dependence on the Majorana mass term M(0), 
which results in p{4>) = exp [— 3(a + 1)0/2]. The third case describes a sterile neutrino with 
a dynamical mass evolving like the GCG with mass varying dark matter coupled to dark 



M(4>) = M 



tanh (3(a + l)f) 
tanh (3(a + l)f ) 
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energy with a equation of state, p DE = —p DE - Since the prescription for decoupling mass 
varying dark matter from the effective GCG leads to sterile neutrino mass given by Eq. ff28|) , 

2a 

the active neutrino mass should be given by = tanh [3(a + 1)0/2] a+1 exp [— 3(a + 1)0]. 
Keeping in mind the connection established by the seesaw mechanism for the third case, thus 

2a 

the mass matrix term should be given by m{4>) = tanh [3(a + 1)0/2] "+ 1 exp [— 3(a + 1)0/2]. 

Considering that the cosmological background evolves according to the GCG model, 
Fig. [2] exhibits the behavior of the dark sector for the three considered cases. The evolving 
neutrino masses lead to evolution of the neutrino energy densities which are shown in Fig [3j 
where active and sterile neutrino energy density ratios, p A /ps, are plotted as function of 
the scale factor a. The sterile neutrino is assumed to behave like dark matter. Notice that 
different thermodynamic regimes are considered at present: NR, relativistic and UR. As 
expected, once the active neutrino reaches the NR regime, the masses and densities are not 
independent and follow the exponential behavior given by p A j p s oc s 2 = exp [3(a + 1)0]. 

The equation of state for active neutrinos in NR regime at present is shown in Fig. HI 
Despite the explicit model dependence of the masses on the scalar field, Figs. |3] and H] show 
similar global features that are typical of the proposed framework. The resulting mass 
dependence on the scale factor for the three considered cases is depicted in Fig. [5j In the 
three cases one sees a sharp increase of the neutrino masses at a recent past (a > 0.2). 

Our analysis shows that the active neutrino mass component grows when coupled to dark 
matter through the seesaw cosmon field within the GCG model. This is consistent with the 
treatment of neutrinos as a test fluid, in the sense that its cosmological equation of motion 
decouples from the dark sector which governs the cosmological evolution. In addition, it 
provides the conditions to analyze instabilities. The possibility of adiabatic instabilities in 
cosmological scenarios was previously pointed out in the context of a mass varying neutrino 
model of dark energy whereas the dynamical dark energy model is obtained through the 
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coupling of a light scalar field to neutrinos, but not to dark matter 
have been extensively discussed in the context of mass varying neutrinos, in which the growth 
of the neutrino's mass and the recent accelerated expansion are linked through the scalar field 
coupling. In the adiabatic regime, these models face catastrophic instabilities on small scales, 
since the sound speed squared of the coupled fluid is negative. Starting with a uniform fluid, 
such instabilities would give rise to exponential growth of small perturbations. The natural 
interpretation of this is that the universe becomes inhomogeneous with neutrino overdensities 
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subject to nonlinear fluctuations which eventually collapse into compact localized regions 

A 

Our approach allows one to perform a stability analysis for NR active neutrinos indepen- 
dently of the dark sector evolution. The first step is obtaining the squared speed of sound, 
c 2 s , for neutrinos. 

Fig.[6]illustrates the results for c\ for DFG of neutrinos {y), the seesaw cosmon dark energy 
and the GCG plus neutrinos. The results are compared with those of a GCG scenario with 
A 3 = 4/5 for a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot line), and 1/8 (dash line). 
The role of neutrinos on the positiveness of c 2 s is measured in terms of those for the GCG 
scenario. One sees that the influence of neutrinos on the positiveness of c 2 s for the GCG is 
not relevant. Actually for the first two cases, NR neutrinos lead to vanishing perturbations 
on the background at present, so that the test fluid approach is quite accurate. 

The results obtained from the analysis of the third case show that for the models where 
the stationary condition (cf. Eq. ( )27l) ) implies that, = —p^,, one obtains c 2 s = — 1 from 
the very start. The effective GCG plus neutrinos is free from this inconsistency, i. e. the 
coupling of the dark energy component with mass varying dark matter allows for removing 
such inconsistency as c 2 s ~ dp^/dp^ > 0. 

VI. CONCLUSIONS 

In this work it is shown that GCG contains some key ingredients in the description of the 
universe dynamics and that it allows for a suitable background for the mass varying and the 
seesaw mechanisms. Furthermore it has been pointed out that the mass varying behavior of 
the dark matter component can be matched with the GCG model. 

The main features of the proposed model can be summarized as follows: i) The cosmolog- 
ical evolution of neutrino energy densities is driven by coupled cosmon-type field equations 
where active and sterile neutrino states are connected through the seesaw mechanism for 
mass generation, ii) Dark matter is, most often, not considered in the MaVaN models. The 
treatment of dark energy and dark matter in the GCG unified scheme naturally offers this 
possibility. Identifying sterile neutrinos as dark matter coupled with dark energy provides 
the conditions to implement such unified picture in the MaVaN formulation. Moreover, 
the constraints imposed by the seesaw mechanism in order to establish a unique analytical 
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connection to the GCG in terms of a real scalar field were found, iii) The dynamics of the 
coupled fluid composed by neutrinos, dark matter and dark energy is driven by one single 
degree of freedom, the scalar field, 0(a). Since the GCG allows for an explicit dependence 
on the scale factor, a, s and £ dependence on the scalar field do imply a dependence on 
the universe's evolution. Due to the connection between the GCG and the seesaw masses, 
the proposed approach actually yields an effective model for MaVaN's coupling to the dark 
sector. Of course, one can interpret the mediation of the scalar field as a dependence on the 
dimensionless scales a/a , s, and £/M. 

The described mechanism leads to a fast increase of the neutrino masses, which results in a 
model dependent vanishing speed of sound at present. The dynamical mass behavior, due to 
the evolution of the seesaw cosmon field, coincides in a subtle way with the GCG dynamics. 
However, this scenario is not unique. Nevertheless, without any additional assumptions, 
it provides an attractive solution for the coincidence problem through the confluence of 
three independent frameworks: the cosmon-like dynamics, the seesaw mechanism for mass 
generation and the GCG model. For active neutrinos, an increase of the mass /i by a factor 
10 6 , corresponds approximately to a decrease of the sterile neutrino mass Ai from the Planck 
scale to 10 13 GeV. Thus, the GCG model modulates the increase of the neutrino mass, which 
acts as a cosmological clock to the present scenario of cosmological accelerated expansion. 

Unfortunately, in general, one cannot provide a sharp criterion for the potential and for 
the mass varying dependence on the scalar field to discriminate between the GCG scenario 
and a unified fluid via a seesaw cosmon field that mimics the GCG. Furthermore, it has 
been shown for specific potentials that, in what concerns stability and cosmic acceleration, 
many results found in the literature can be recovered. Actually, new scenarios featuring 
other mass dependencies on the seesaw cosmon field can also be considered. 

To conclude, the proposed unified scheme, despite similarities with some quintessence 
models, is an ambitious and encompassing scheme where mass varying particles can be 
related with the stability and the cosmic acceleration of the universe. 
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Appendix: The mass varying mechanism for a DFG 



Aiming to understand how the mass varying behavior takes place, it is assumed that for 

m 

T 



a particle statistical distribution f(q), where q = T being the background temperature 



at present, in a flat FRW cosmo logical scenario. For a generic dynamical mass defined by 
M (4>), the corresponding particle density, energy density and pressure are expressed as 



11(a) 
p(a,(f>) 
p(a,4>) 



± 



T'4 
n 



T 4 
1 o 

37r 2 a 4 



dq q 2 f(q), 



dq q 2 [ q 2 + 



M 2 (<j>) a 2 V /2 



rp 2 
X 



dq g 4 [q 2 + 



M 2 Q) a 2 

^ 



/(?), 

-1/2 



(29) 



/(<?), 



where one has introduced the sub-index for denoting present-day values, with a = 1. 

In the limit where T tends to 0, the Fermi distribution f(q) becomes a step function that 
yields an elementary integral with the upper limit equal to the Fermi momentum qp. The 
energy density and pressure of a DFG can be expressed in terms of elementary functions of 
the scale factor, (3 = (3(a) = T q F /a and M = M(</>(a)), 



n(a) = 


3vr 2 ^ ' 
1 r 


p(a) = 




(8tt 2 ) I 
1 


p(a) = 


(8tt 2 ) 



/?(2/3 2 + M 2 ) + M 2 - M 4 arc sinh (/3/M) 



/3(-/3 2 - M 2 )v / ^ 7 TM 2 + M 4 arc sinh (/3/M) 
3 



Simple mathematical manipulations allow one to show that 

dp(a) 



(30) 



Tl(a) 



and 



M(o) 



dn(a) 
dM(a) 



(p{a)+p(a)), 



(p(a) - 3p(o)), 



from which one can obtain the energy-momentum conservation equation 

P + 3tf( P + p)-0^ = O. 



(31) 



(32) 



(33) 



The coupling between relic particles and the scalar field as described by Eq. ( |32l) is relevant 
only in the NR regime since the strength of the coupling is suppressed by the pressure of the 
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relativistic (T(a) = T /a » M(<j>(a))) particles, so that matter and scalar field fluids tend to 
decouple and evolve adiabatically. 
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FIG. 1: The correspondence between the energy density components and the "modified" seesaw 
scale parameter 7s 2 (first plot), and the equation of state Wi = Pi/pi- Sterile and active neutrinos 
are assumed to behave like a DFG in NR (77 = 0.1), relativistic (77 = 1) and UR (77 = 10) regimes. 
For the NR limit of a DFG, on has Pa/Ps ~ 7 _3 s 2 (straight solid line in the plot). 
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FIG. 2: Energy densities £1 = p(a)/p in terms of the scale factor, a (with a = 1), for the 
components of a unified background fluid (effective GCG): mass varying dark matter (DM) 
p s and cosmon-like dark energy (DE). Three cases are considered for the dynamical neutrino 
masses: p(4>) = 0exp [— 3(a + l)t/>/2] (Case 01), p((j)) = exp [— 3(a + l)</>/2] (Case 02) and p((j>) = 
tanh [3(a + l)<fi/2] "+ 1 exp [— 3(a + l)4>} (Case 03). The dark matter densities are computed for M.(4>) 
in correspondence to the GCG scenario through Eq. (|23p . with A s = 4/5 and a = 1 (solid line), 
1/2 (dot line), 1/4 (dash dot line), and 1/8 (dash line). 
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FIG. 3: Active and sterile neutrino energy density rate p A j ' ps as function of the scale factor, a 
(with a = 1). Both sterile and active components are assumed to behave like a DFG. The results 
are obtained for different GCG models with A s = 4/5 and parameters a = 1 (solid line), 1/2 
(dot line), 1/4 (dash dot line), and 1/8 (dash line). Three cases are considered for the neutrino 
masses (Case 01: red lines, Case 02: blue lines, Case 03: green lines) in agreement with Fig. [2] for 
nonrelativistic (NR), relativistic and ultrarelativistic (UR) neutrinos at present. 
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FIG. 4: Neutrino equation of state w = Pa/ Pa for the GCG cosmological model with neutrinos 
behaving like a DFG in the NR regime (77 = 0.1 at present). The GCG parameters and the neutrino 
mass dependence on 4> are in correspondence with the previous figures. 
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FIG. 5: Growing neutrino mass (//(a)//x ) i n dependence on the scale factor, a. Masses are nor- 
malized in terms of fj, at present (o = o = 1). The GCG parameters and the neutrino mass 
dependence on <fi are in correspondence with the ones in the previous figures. The GCG scenarios 
where A 3 = 4/5 and a = 1, 1/2, 1/4, 1/8 are considered. 
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FIG. 6: Squared speed of sound c 2 s = dp /dp in dependence on the scale factor, a (with a = 1), for 
neutrinos (magenta lines), for the cosmon-like dark energy component - DE (black lines) and for the 
corresponding GCG scenario (blue lines). The GCG parameters and the neutrino mass dependence 
on (j> are in correspondence with the previous figures. Neutrinos are assumed to behave like a DFG 
in the NR regime (rj = 0.1 at present). 
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